About the fastest growth of Order Parameter in Models of Percolation 
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Can there be a 'Litmus test' for determining the nature of transition in models of percolation? 
In this paper we argue that the answer is in the affirmative. All one needs to do is to measure the 
'growth exponent' x of the largest component at the percolation threshold; \ < 1 or \ — 1 determines 
if the transition is continuous or discontinuous. We show that a related exponent rj = 1 — \ which 
describes how the average maximal jump sizes in the Order Parameter decays on increasing the 
system size, is the single exponent that describes the finite-size scaling of a number of distributions 
related to the fastest growth of the Order Parameter in these problems. Excellent quality scaling 
analysis are presented for the two single peak distributions corresponding to the Order Parameters 
at the two ends of the maximal jump, the bimodal distribution constructed by interpolation of these 
distributions and for the distribution of the maximal jump in the Order Parameter. 

PACS numbers: 64.60.ah 64.60.De 64.60.aq 89.75.Hc 



Macroscopic properties change when a thermodynamic 
system undergoes a phase transition from one phase to 
another. The transition is characterized by looking at 
how a certain observable, called the Order Parameter 
(OP), varies before, at and beyond the critical point of 
a transition while a suitable control variable is contin- 
uously tuned. The OP exhibits a very rapid change in 
magnitude starting from the critical point. Depending 
on if this rapid change occurs discontinuously or contin- 
uously, the transition is termed as discontinuous (or first 
order) or continuous (or second order) transition [1]. 

Recently it has been seen that such straight forward 
classification may not be easy in some models of percola- 
tion phenomena @,I3. Originally the percolation model 
was introduced by Broadbent and Hammersley in 1957 
0] to better understand the mechanism of charcoal based 
gas masks and later this model had been used exten- 
sively to study the order / disorder transitions in vari- 
ous systems. Long range correlation in terms of global 
connectivity appears as the density of nodes / links is 
increased beyond certain critical point, called the per- 
colation threshold. Over several decades many different 
variants of percolation model had all exhibited only con- 
tinuous transitions. This was the situation until recently 
when Achlioptas et. al. argued that a local competition 
between a pair of vacant edges for being occupied indeed 
leads to an abrupt jump in OP. Hence they coined the 
name 'Explosive Percolation' (EP) [|| to emphasize the 
abruptness of the transition. 

The original model of EP had been defined by a slight 
modification of Random Graphs [6] and we refer to it as 
Achlioptas process (AP) in the following. One starts with 
N nodes (each node is a component of size unity) and 
no links. At an arbitrary intermediate stage a pair of va- 
cant edges are randomly selected between nodes and 
(k, I) which belong to the components of sizes S{, Sj and 
Sfc, si respectively. There is a competition: if SiSj < SfcS; 
then link is occupied, otherwise (k,l) is occupied; 

however when the products are equal, one of the two 



edges is selected randomly. This bias towards connect- 
ing small components delays the growth of the largest 
component and consequently the percolation threshold 
is pushed up to p c = 0.888449(2) Q compared to 1/2 for 
Random Graphs. In general the preferential link occupa- 
tion rule slows down the growth of largest component and 
accelerates smaller components to grow faster. Further 
EP has been studied on the Square Lattice [8j,|9j, Random 

Graph [ltj, scale- free networks [nun, with 

a Hamilto- 



nian formulation [14|, in a cluster aggregation process 



15 1 and also in human protein homology network [l6j |. 
real-world networks [17] etc. 

Subsequently situation started changing and a num- 
ber of research publications suggested that EP transition 
is actually continuous. The first claim came from da 
Costa et. al. who considered a slightly modified ver- 



sion of AP [18|| . For each of the randomly selected pair 
of vacant edges one first picks up the node with smaller 
component size and then connects them. It was argued 
that though this model is more biased than AP it has 
a continuous transition and therefore the original AP 
must then have a continuous transition. Later it was 
shown numerically that the magnitude of the maximal 
jump in OP has a zero measure in the thermodynamic 
limit [l|| [2(| . Very recently Riordan and Warnke have 
rigorously proved that all models using Achlioptas type 
processes are indeed continuous in the asymptotic limit 

A distinction between two types of EP models can be 
made: (i) Local rule: As in Random Graphs these mod- 
els pick up in general a few nodes randomly with uniform 
probability and then preferentially occupy a vacant edge 
within this subset of nodes to introduce the bias towards 
connecting smaller components, (ii) Global rule: The 
entire set of vacant edges are assigned non- uniform prob- 
abilities which are already biased towards smaller com- 
ponents and out of them one vacant edge is randomly 
selected for occupation. In such a model a vacant edge 
is selected with a probability proportional to {stSj)^ and 
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FIG. 1: (Color online) The AP data of a single run. (a) 
The jump sizes AC(p,N) in the Order Parameter have been 
plotted (blue symbols) with link density p for the graph of 
size N = 10 6 . The subset of red symbols corresponds to those 
jumps where existing largest components increased their sizes, 
(b) The average jump size {AC(p, N)} plotted with deviation 
p —p c (N) from the percolation threshold for N = 10 5 . 
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FIG. 2: (Color online) For AP the finite-size scaling of av- 
erage size of the jump in OP for N = 10 4 (blue), 10 s (red) 
and 10 6 (green), (a) Sub-critical region. The average jump 
size (AC<(p, N)) is plotted with deviation \p — p c (N)\. The 
average slopes give q< = 1.18(5). (b) Super-critical region. 
The average jump size (AC>(p, N)) is plotted with deviation 
p — p c (N). The average slopes give a> = 1.02(5). 



it is argued that for all values of £ < £ c the transition is 
discontinuous [19j. 

The distribution of the Order Parameter at the per- 
colation threshold has a bimodal distribution where the 
depth of valley between the two peaks vanishes as N — > 
oo which signifies a discontinuous transition 

0,113, El. 

Extensive scaling analysis of these distributions in a set 
of four Local rule models showed that all of them have 
continuous transitions Q- The positions of peaks scale 
with N with different exponents which indicate different 
growth mechanisms for the sub and super-critical perco- 
lation regimes Q- 

In all percolation models we have studied, we start 
from the empty graphs with N nodes and no edges. Links 
are then dropped one by one following the specific rules of 
the model and sizes of different components of the graph 
start growing. The link density is defined as p = n/N 
where n is the number of links in the graph. If s max 
is the size of the largest component of the graph then 
the Order Parameter is defined as C(N) = s max /N. To 
characterize the percolation transition we study the per- 
colating system when the OP undergoes its fastest rate of 
growth. As links are being dropped at unit rate we keep 
track of s max and calculate how much s max increases due 
to each link addition. In a certain run, let the maximal 



jump in Smax occur due to addition of n m ax-th link. We 
then define the percolation threshold of this specific run 
as Pmax(N) = n max /N. Consequently the percolation 
threshold p c {N) of a graph of size N has been determined 
by the average of p max (N) values over a large number of 
runs, i.e., 



Pc{N) = (p ■max 



(i) 



It is first noted that there is hardly any unique largest 
component below the percolation threshold. More specif- 
ically addition of a single link may merge two smaller 
components whose combined size may be larger than 
Smax- To get a first hand idea we plot in Fig. 1(a) 
using blue symbols the successive jump heights in OP 
i.e., AC(N) vs. p in a single run for AP with N = 10 6 . 
Further we re-plot a subset of these points by red sym- 
bols for whom the jumps in s m ax occurred to the al- 
ready existing largest component. This plot shows that 
the percolation threshold p max (N) clearly distinguishes 
two regimes. In the sub-critical regime p < p m ax(N) 
there is few red points, signifying a unique largest com- 
ponent of the graph does not exist in this regime of link 
density. On the other hand in the super-critical regime 
P > Pmax{N) the jumps occur much more frequently than 
in the sub-critical regime since the largest component is 
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FIG. 3: (Color online) The average values of (C L (N)) (filled 
circles) and {Ch(N)) (empty circles) are plotted with TV for 
AP (blue) and for da Costa model (red). From slopes the 
growth exponent \ °f the largest component are estimated as 
0.9355(5) for AP and 0.9553(5) for da Costa model. 

already so big that almost any link addition enhances 
its size. There cannot be a blue point in this regime 
since that would imply the appearance of another largest 
component, which will eventually merge with present one 
with an even larger jump in C which is not possible since 
the maximal jump in C has already occurred. 

In Fig. 1(b) we plot the average jump size (AC(p, TV)) 
with the deviation p—p c (TV) from the percolation thresh- 
old. The variation of AC(p,N) is different in the sub- 
critical regime and in the super-critical regime. This 
looks similar to the well known A-transition of the diver- 
gence of specific heat in critical phenomena characterized 
by the same value of exponent but with different ampli- 
tudes for the sub and super-critical regimes. In compar- 
ison we get different exponents as reported below for the 
two sides of the percolation threshold. 

In Fig. 2 this analysis has been done in more de- 
tail and a finite-size scaling analysis of (AC(p, TV)) is 
shown. In the sub-critical regime the average jump 
size (AC<(p, TV)) have been calculated as a function of 
\p — p c (N)\. For all three system sizes the plots of binned 
data on a log — log scale exhibit similar behavior, i.e., an 
initial horizontal part for very small values of \p — p c (N)\ 
followed by a linear downward regime, indicating a power 
law decay like: 

<AC<(p,TV)) ~ \ P - Pc (N)\~ a < . (2) 

The directly measured values of a< are 1.133, 1.157 and 
1.164 respectively for the three system sizes which ap- 
proach to 1.18(5) for large TV. In Fig. 2(a) a finite-size 
scaling analysis of this data has been done using the form: 

<AC<(p,TV))TV - 25 ~e<[|p-p c (iV)|/V - 70 ]. ( 3 ) 

For the super-critical regime the average jump size 
(AC>(p,JV)) in OP has been plotted with p - p c (N) in 
Fig. 2(b). A similar variation of initial flat, followed by 
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FIG. 4: (Color online) (a) The average maximal jump 
(ACmax(N)} in OP plotted on a log — log scale with graph 
size TV. The slopes are 7j(AP)f» 0.06475 (red) and 7?(da 
Costa)~ 0.04495 (blue), (b) The same sets of data are plotted 
against N~ v on a lin-lin scale. Continuous lines are straight 
line fits of the data which are then extrapolated to TV — s> oo to 
meet (AC roTC (TV)) axis at {AC max (oo)) « -0.000516 for AP 
(red) and « -0.000256 for da Costa model (blue). 

a power law decay {AC > (p, TV)) ~ (p — p c (N))~ a> has 
been found. The directly measured values of a > are 0.79, 
0.89 and 0.92 respectively which approach to 1.02(5) for 
large TV. A similar data collapse was possible using the 
following scaling equation: 

<AC>(p,TV))TV - 27 - Pc (N))N - 77 }. (4) 

It seems likely that growth rules of few leading com- 
ponents are similar at the percolation threshold. For ex- 
ample the largest component, second largest component 
etc. may grow in the same way since they grow 'indepen- 
dently' but in very similar environments. We define the 
'growth exponent' x f° r the size of the largest compo- 
nent (and next few leading components) at the percola- 
tion threshold as (s max {j) c (N), TV)) ~ TV X or in terms of 
OP (C(p c (TV), TV)) - TV - * 7 with r\ = 1 - x- Since for reg- 
ular d-dimensional lattices TV = L d , percolation on these 
lattices has the growth exponent x — D/d where D is 
the fractal dimension of the 'infinite' incipient cluster at 
the percolation threshold. 

While the size of the largest component monotonically 
increases with link density size of the second largest com- 
ponent gradually increases, reaches a maximum and then 
suddenly falls to a small value [23j |. The maximal jump 
in the OP occurs precisely at this point when the largest 



4 




0.1 0.2 0.3 0.4 0.5 0.6 0.00 0.25 0.50 0.75 1.00 1.25 0.25 0.50 0.75 1.00 1.25 1.50 




FIG. 5: (Color online) For Achlioptas Process on graphs of sizes N = 10 3 (black), 10 4 (red) and 10 5 (blue). Finite -size scaling 
analysis are displayed for different probability distributions associated with the Order Parameter C using rj(AP)= 0.0645. They 
are: (a) Maximal jump sizes ACmax(N). (b) The lower end value Cl{N) of the maximal jump in C. (c) The higher end value 
Ch(N) of the maximal jump in C. (d) The interpolating function C W (N) of Cl(N) and Ch(N) with probabilities r and 1 — r. 
The right-top curve is for r=0.1 and increased to 0.9 at the interval of 0.1 using only one system size N = 10 3 . (e) The scaling 
of P(C W (N)) for r = 0.3 (right-top), 0.5 and 0.7. (f) The ratio of heights of the right peak hn(r) and left peak hh{r) of the 
bimodal distributions in 5(d) plotted with probability r. 



component merges with the largest (in the entire his- 
tory of the system) second largest component. Since 
both of them grow in a similar way it is expected that 
the height of the maximal jump in the largest compo- 
nent should also grow as N x . Let us assume that dur- 
ing the fastest growth OP jumps from Cl{N) to Cr{N) 
so that AC max (N) = C H {N) - C L (N). In Fig. 3 we 
plot on a log — log scale the average values of {Cl(N)) 
and (Ch{N)) with N for AP. Linear least square fits ex- 
hibit that both plots are highly straight, no systematic 
variation of slopes could be detected and they are very 
closely parallel to each other with slopes -0.0641(5) and 
-0.0645(5) giving an average of ^(AP)=0. 0643(5) and the 
growth exponent x(AP)= 0.9357(5). Similar data for the 
da Costa model have also been plotted in Fig. 3 with red 
color and the estimates for the slopes are -0.0443(5) and 
-0.0445(5) corresponding to the growth exponent of x(da 
Costa)= 0.9556(5). 

In addition, size of the maximal jump has also been es- 
timated directly and its average value for different system 
sizes are then extrapolated as: 

(AC max (N)) = (AC max (oo)) + AN'* (5) 

as N — > oo and if (AC ma :r(oc)) turns out to be zero or 
non-zero it signifies a continuous / discontinuous change 
in OP. In Fig. 4(a) we plot (AC max (N)) vs. N for AP us- 



ing a log — log scale. A straight line fits to the data very 
accurately except for very small system sizes and gives 
r;(AP)= 0.06475. In Fig. 4(b) we use this value of r\ and 
plot (AC max {N)) vs. N~ v on a lin - lin scale. Here again 
the least square fit of a straight line works nicely, giving 
(AC max (oo)) = —0.000516. This analysis suggests that 
indeed the form of Eqn. (4) is very likely to be valid and 
since the constant (AC max (oo)} is nearly equal to zero, 
it gives a strong indication that the transition of AP is 
actually continuous. Similarly for da Costa model we get 
?7(da Costa)= 0.04495 and (AC max (oo)) = -0.000256. 
We conclude r/(AP)= 0.0645(5), x(AP)= 0.9355(5); rj(d& 
Costa)= 0.0447(5), x(da Costa)= 0.9553(5) and use 
these values in the following analysis. 

Fig. 5 shows the plots for the finite-size scaling anal- 
ysis of a number of probability distributions associated 
with the maximal jump of OP in the Achlioptas Process. 
They are for: (a) the maximal jump size AC max (N); (b) 
the lower end value Cl{N) of the maximal jump; (c) the 
higher end value Ch(N) of the maximal jump. While the 
distribution P(Cl(N)) appears to be almost symmetric 
the distribution P(Ch(N)) is found to be quite asym- 
metric. For the scaling analysis three system sizes are 
used in each case and the data collapse are found to be 
excellent. In each case the x-axis is scaled by N^ Ap ) 
and the y-axis is scaled by N~ v ^ Aps> . We further define 
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FIG. 6: (Color online) Distribution of the Order Parame- 
ter C(p c ,N) right at the percolation threshold p c (N). (a) 
P{C(p c ,N)) plotted with C(p c ,N) for N = 10000 (black), 
32000 (red) and 100000 (blue). Finite-size scaling of the right 
peak in (b) and that of the left peak in (c). 



a single value C w of the Order Parameter at the maxi- 
mal jump as C w — Cl with probability r and C w = Ch 
with probability 1 — r. One can define the interpolat- 
ing function with respect to the common variable C as: 
P r (C) = rP L {C) + (l-r)P H (C). In Fig. 5(d) we show this 
interpolating function P(C W (N)) having double humps 
with unequal heights and widths for 9 values of r from 0.1 
to 0.9 at an interval of 0.1 for only one size 7V=1000. In 
Fig. 5(e) we plot the scaling analysis of P(C W (N)) for few 
representative values of r = 0.3, 0.5 and 0.7 which also 
scale with N r, ^ AP \ In both Fig. 5(d) and 5(e) the com- 
mon point of intersection where all curves pass through 
has coordinates (0.855, 1.14). Finally in Fig. 5(f) we plot 
the ratio hn(r)/hL(r) of the heights of right peak and 
left peak as a function of r. Therefore it turns out that 
only one exponent ?y(AP) which determines the growth 
exponent x(AP)= 1 — 77 (AP) characterizes various dis- 
tributions related to the fastest growth of the Order Pa- 
rameter in AP. In addition it is also verified that for da 




FIG. 7: (Color online) Plot of p c (N) - p c (oo) vs. N on the 
log — log scale for AP (black), da Costa model (red) and Ran- 
dom Graphs (blue), (a) Apart from few small size graphs the 
data points fit nicely to straight lines giving values of the v 
exponent in Eqn. (6) as 1.50, 1.26 and 2.98 respectively, (b) 
Similar data for the AP (empty circles) and ordinary bond 
percolation (filled circles) on square lattice giving v — 1.98 
and 2.67 respectively. 



Costa model all these probability distributions are very 
much similar and they scale equally well with the corre- 
sponding value of 7y(da Costa). 

Bimodal distributions have been observed for the prob- 
ability of the Order Parameter at the percolation thresh- 
old |7J, [l3|, 122). However this distribution depends on 
the precise definition of the percolation threshold. In [?| 



the percolation threshold p' c (N) has been determined in 
such a way so that heights of the two peaks are equal 
and this is kept fixed for all runs. In comparison in our 
case different runs have different percolation thresholds 
Pmax(N) since the maximal jumps AC max {N) in OP oc- 
cur at these specific link densities. These are two different 
statistical ensembles though they approach the same per- 
colation point p c (oo) as N — > 00. Evidently in a single 
run Pmax(N) may be either smaller or larger than p' c (N). 
If p' c (N) < p m ax{N) then its corresponding OP C < Cl 
and when p' c (N) > p ma x(N) then its C > Ch- For this 
reason with a fixed run-independent value of p' c (N) there 
are two regions where C occur more frequently than its 
other values resulting a bimodal distribution of OP. The 
peaks of the bimodal distribution occurring at C± with 
C+ > C_ scale as C± ~ N~ v± 0. We observe that our 
exponent rj lies in between them i.e., rj + < r/ < r;_ for 
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FIG. 9: (Color online) For the Global rule EP model on 
square lattice the values of (Cl(N)) (square), (Ch(N)) (trian- 
gle) and {ACmax(N)) (circle) are plotted for £ = -1 (black), 
-2(red) and -3(blue). 
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FIG. 8: (Color online) For AP the probability distributions 
P(Pmax(N)) vs. Pmax(N) — p c (N) for three system sizes N 
! I! ' (black), 10 4 (red) and 10 5 (blue) in (a) and their finite-size 
scaling in (b). 

both AP and da Costa model (Table I). 

Further we calculate the distribution of the Order Pa- 
rameter at p c (N) defined in Eqn. (1) and fixed for all 
runs. A sample size of 10 7 has been used for three differ- 
ent system sizes N = 10000, 32000 and 100000. These are 
bimodal distributions P(C Pc (N)) with unequal heights 
and widths shown in Fig. 6(a). We see positions of both 
peaks C±(p c , N) decreases as N increases. The finite-size 
scaling of the right and left peaks are displayed in Fig. 
6(b) and Fig. 6(c). Assuming the same scaling forms 
C±(p c , N) ~ N~ r > ± we obtain -q+ = 0.042 and r)- = 0.31. 
Here also the same inequality r/ + < rj < r/_ holds good. 

Next we make a systematic study how the percola- 
tion threshold p c (N) approaches to its asymptotic value 
p c {oo). We cannot talk about a correlation length expo- 
nent in this case since graphs are not embedded in Eu- 
clidean space, yet we assume a leading order correction 
for the finite size systems associated with an exponent v 
(not the correlation length exponent) as: 

p c (N) = p c (oo) + BN~ 1/u . (6) 

In Fig. 7(a) we plot p c {N) — p c {oo) vs. iVona log — log 
scale using p c (oo) — 0.888446 (note that this differs 3 
in the last digit from the value quoted in |7j) for AP, 
0.923207508 for the da Costa model [l8] and 1/2 for 
the Random Graphs Q. On such a scale each plot is 
expected to be linear if the corresponding p c {N) val- 
ues obey the Eqn. (6). It is observed that apart from 



very small graphs the data points do fit to straight lines. 
The slopes in these regions give v values as 1.50, 1.26 
and 2.98 for AP, da Costa model and Random Graphs 
which we conjecture may be 3/2, 5/4 and 3 exactly. We 
further plot similar data of AP on the square lattice in 
Fig. 7(b) using p c (oo) = 0.526565 Q. The curve is a 
very nice straight line with slope -0.505 corresponding to 
v ss 1.98 which leads us to conjecture vap = 2 exactly 
for square lattice. With this result it seems to be bet- 
ter to fit p c {N) data for square lattice of size L (with 
N = L 2 ) as p c {N) = p c {oo) + ciA^ 1 / 2 + c^" 1 [H. 
Using this quadratic fitting form and with eight data 
points from L = 32 to 4096 we get p c (oo) = 0.526575. 
The extrapolated value systematically came down slowly 
when we discarded the data for small lattices. Finally 
for the last three points L = 1024, 2048 and 4096 we 
get Pc(oo) = 0.5265639. In addition we repeat this cal- 
culation for ordinary bond percolation on square lattice 
and using p c (oo) = 1/2 we get a good straight line with 
v w 2.67 (Fig. 7(b)). 

Next we calculate for AP the probability density dis- 
tribution of the percolation thresholds P(p m ax{N)) and 
plot them with p ma x{N) — p c {N) in Fig. 8(a) again for 
the same three different system sizes and with 10 7 con- 
figurations for each system size. All three plots fit very 
closely to Gaussian distributions with different parame- 
ter values for different N. Fig. 8(b) shows the finite-size 
scaling of these data and the best data collapse is ob- 
tained corresponding to the following form: 

P{p max {N)) ~ N 9 T([p max {N) ~p c {N)]N e ). (7) 

A best value of 9 = 0.51 is obtained which is very much 
consistent with the conjecture of 9 = 1/2 in Q. 

Finally we calculate the growth exponent x again for 
the Global rule EP model [19(. It was argued that this 
model exhibits discontinuous transitions for all values of 
the biasing parameter ( < on the square lattice. We 



7 





AP 


da Costa 


AP in 2d 


BP 


RG 


Vc 


0.888449(2) [7] 


0.923207508 [18] 


0.526565 [9] 


1/2 


1/2 




0.888446 




0.5265639 






n+ 


0.0402(15) [7] 


0.0255(80) [7] 


0.018(2) [7J 






V 


0.0645(5) 


0.0447(5) 


0.0217(5) 


5/96 


1/3 


V- 


0.270(7) [7] 


0.300(5) [7] 


0.078(7) [7] 






X 


0.9355(5) 


0.9553 


0.9783(5) 


91/96 


2/3 


V 


1.50 


1.26 


1.98 


8/3 


2.98 



TABLE I: Values of different exponents available in the liter- 
ature as well as measured in this work. Some known results 
of Bond Percolation (BP) in 2d and for Random Graphs are 
also included for comparison. The p c , r/+ and r\- values are 
taken form Q. 



therefore calculate (C L (N)), (C H (N)) and (AC max (N)) 
for ( = —1, —2 and -3. In Fig. 9 we plot them with N on 
a lin-log scale and see that curves are simply horizontal 
straight lines indicating N independence. This implies 
that the growth exponent \ = 1 and r\ = for the Global 
rule EP model [H. Therefore in general (AC max (N)) = 
constant implies that for these C values the transitions 
are discontinuous. 

We summarize that by looking at the numerical value 
of the exponent x describing the growth of the largest 
component at the percolation thresholdwith with graph 
size N one should be able to understand if the model 
exhibits a continuous or discontinuous transition. The 
growth exponent is defined as (s max ) ~ N x has been 
quite well known for graphs in the literature but perhaps 
without a name. This exponent is similar to the frac- 
tal dimension when the graph is embedded in Euclidean 
space i.e., X < 1 implies a 'fractal' with fractal dimension 
being less than the embedding Euclidean space dimension 
and % = 1 implies a 'compact' giant component whose 
size grows proportional to the graph size N. We justify 
this as follows. The maximal jump in the Order Parame- 
ter takes place only when the largest component merges 
with the largest second largest component and sizes of 
both grow as N x at the percolation threshold. For a 
continuous transition the maximal jump in OP must have 
a zero measure with respect to the graph size N in the 
asymptotic limit and therefore % < 1. On the other hand 
for a discontinuous transition, the largest jump in OP 
which is the size of the largest second largest component 
must be proportional to the graph size N and therefore 
X = 1. We conclude that any percolation model which 
has x < 1 / X = 1 must exhibit continuous / discontinu- 
ous transition and the reverse is also true. In addition we 
show that a number of distributions related to the largest 
jump of OP at the percolation threshold are described by 
a related exponent r\ = 1 — x- A comparison of all related 
exponents have been done in Table 1. 

In recent times it is being shown that asymptoticall y al l 
Local rule models have continuous transitions 0, flii - E^ . 
In comparison Global rule models exhibit discontinuous 



transitions. We conjecture that for all Local rule mod- 
els with x < 1 transitions are continuous, whereas all 
Global rule models have x — 1 an d their transitions are 
discontinuous. 

Finally we would like to make the following comment. 
Melting of ice is a well known example of first order tran- 
sition where the density of ice changes from 0.92 gm/ml 
to 1.00 gm/ml of water. On the other hand on a graph 
with Avogadro number of nodes (N ~ 10 23 ) the Achliop- 
tas process has the largest jump in the Order Parmeter 
as AC max {N) = N-^ AP ^ » 0.03 with r/{AP) = 0.0645. 
Therefore we believe that 'practically' the transition in 
Achlioptas process may very well be considered as the 
discontinuous transition. 

I am very much indebted to Deepak Dhar, Peter Grass- 
berger, Maya Paczuski and Robert M. Ziff for many help- 
ful suggestions. I convey my sincere thanks to Arnab 
Chatterjee, Raissa D'Souza and Janos Kertesz for the 
critical reading of the manuscript. 
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